Hyperpolarizability effects in a Sr optical lattice clock 
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We report the observation of the higher order frequency shift due to the trapping field in a 87 Sr 
optical lattice clock. We show that at the magic wavelength of the lattice, where the first order term 
cancels, the higher order shift will not constitute a limitation to the fractional accuracy of the clock 
at a level of 10 -18 . This result is achieved by operating the clock at very high trapping intensity 
up to 400kW/cm 2 and by a specific study of the effect of the two two- photon transitions near the 
magic wavelength. 
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The recent proposal and preliminary realizations of op- 
tical lattice clocks open a promising route towards fre- 
quencystandards with a fractional accuracy better than 
A large number of atoms are confined 
in micro-traps formed by the interference pattern of laser 
beams which in principle allows both the high signal to 
noise ratio of optical clocks with neutral atoms [5| and 
the cancellation of motional effects of trapped ion de- 
vices 0, • I n contrast to the ion case an optical 
lattice clock requires high trapping fields. The evaluation 
of their effects on the clock transition is a major concern 
and is the subject of this letter. For a Sr optical lattice 
clock, the typical requirement in terms of trapping depth 
is about 10 E r with E r the recoil energy associated to 
the absorption of a lattice photon ^Jj . The correspond- 
ing frequency shift of both clock states amounts to 36 kHz 
at 800 nm, while a relative accuracy goal of 10 -18 implies 
a control of the differential shift at the 0.5 mHz level, or 
10~ 8 in fractional units. 

The frequency of the clock transition in a laser trap 
of depth Uq is shifted with respect to the unperturbed 
frequency vq according to 
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with v\ and proportional to the (dynamic) polariz- 
ability and hyperpolarizability difference between both 
states of the clock transition 1] . By principle of the op- 
tical lattice clock, vi cancels when the laser which forms 
the lattice is tuned to the "magic wavelength" A m . Al- 
though this remains to be demonstrated experimentally, 
a control of this first order term to better than 1 mHz 
seems achievable 0] . 

The higher order term is a priori more problematic 
with no expected cancellation. A theoretical calculation 
of the effect is reported in Ref. 0] predicting a frequency 
shift of —2 fj,Hz/Er for a linear polarization of the lattice. 
The calculation however was performed at the theoretical 
magic wavelength of 800 nm. The actual value [13, A m = 
813.428(1) nm (see 0,0] aim below), lies near two two- 
photon resonances which may considerably enhance the 
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FIG. 1: Energy levels of Sr relevant to this paper. The clock 
transition at 698 nm couples the two lowest energy states of 
the atom. 



effect and impede the realization of an accurate clock. 
The first one couples 5s5p 3 Po to 5s7p 1 P\ (Fig.QJ and 
is at a wavelength of 813.36 nm, or equivalently 30 GHz 
away from the magic wavelength. Although this J — 
— > J = 1 two-photon transition is forbidden to leading 
order for two photons of identical frequencies ^lj , it is so 
close to the magic wavelength that it has to be a priori 
considered. The second one resonantly couples 5s5p 3 Po 
to 5s4/ 3 F 1 2 at 818.57nm and is fully allowed. 

We report here an experimental study of higher order 
effects in a 87 Sr optical lattice clock operating at a very 
high trapping depth up to 1400 E r and for a linear polar- 
ization of the lattice. This depth is about a factor of ten 
higher than in the other reported systems 0,0] which en- 
hances the sensitivity to higher-order frequency shifts by 
two orders of magnitude. At the maximum depth and for 
a detuning of the lattice laser as small as 0.5 GHz from 
the 3 Po — > 1 P\ transition, higher order frequency shifts 
are below a few hertz. On the other hand, when the lat- 
tice laser is tuned to within a few hundred MHz of the 
3 Pq — > 3 F2 transition, non-linear frequency shifts of sev- 
eral kHz are observed. These results, together with mea- 
surements at the magic wavelength, give a stringent limit 
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in optimal operating conditions. We show that higher 
order frequency shifts will not limit the accuracy of a Sr 
optical lattice clock down to the 10 -18 level. 

The high trapping depth is reached thanks to an en- 
hancement Fabry-Perot cavity around the ID vertical lat- 
tice. The circulating power reaches 16 W for a 650 mW 
input at 813 nm. The mode has a waist radius of 90 /im 
corresponding to a maximum Uq = 1400 E r and an ax- 
ial (resp. radial) oscillation frequency of 260 kHz (resp. 
540 Hz). An intracavity dichroic mirror separates the lat- 
tice light from that which probes the clock transition. It 
induces negligible intracavity loss at 813 nm while ensur- 
ing that less than 10 -3 of the incident power at 698 nm 
is reflected back to the atoms by the cavity mirror. The 
lattice polarization is interferometrically filtered by an 
intra-cavity quarter wave plate which lifts the degener- 
acy between the linear polarization states parallel to the 
eigen axis of the plate by half a free spectral range. The 
resulting polarization of the trapping light is linear to 
better than 10~ 3 . Its wavelength Al is controlled by 
means of a wave meter with an accuracy of 10 -3 nm. 

Atoms are loaded into the optical lattice from the 
magneto-optical trap (MOT) described in . Through- 
out the loading cycle the lattice is overlapped with the 
MOT. Cold atoms at the center of the trap are selec- 
tively optically pumped to the metastable 3 Po,2 states 
by means of two " drain" lasers of waist radius 50 /im 
that are aligned to the lattice. They are tuned to the 
- 3 Pi and 3 Pi - 3 Si transitions at 689 and 688 nm 
respectively (Fig.[TJ. Atoms in the metastable states re- 
main trapped provided their energy is smaller than the 
200 /iK lattice depth. This leads to a continuous loading 
of the lattice at a rate of about 10 5 atoms/s. After half 
a second of loading time the MOT and drain lasers are 
switched off, and the atoms are repumped back to the 
ground state using two lasers tuned to the 3 Po,2 — 3 Si 
transitions at 679 and 707 nm. They are then cooled in 
the lattice to ~ 10 /iK in 50 ms with the narrow 1 So — 3 Pi 
intercombination line at 689 nm [l3j | . 

Following this preparation stage we probe the 1 5o — 
3 Po clock transition at 698 nm. The frequency of the 
probe laser is referenced to an ultra-stable cavity as de- 
scribed in Ref . [3| . The probe beam is aligned parallel 
to the lattice and has a waist radius of 200 //m. Its polar- 
ization is linear and parallel to the lattice polarization. 
After the probe pulse, the transition probability to 3 Po 
is measured in a few ms by laser induced fluorescence. A 
first detection pulse at 461 nm gives the number of atoms 
remaining in the ground state and ejects these atoms 
from the trap. Atoms in 3 Po are then repumped to the 
ground state and detected similarly. This method allows 
a transition probability measurement which is insensitive 
to the atom number fluctuations. The resonance is shown 
in Fig.EI The central feature (the carrier) is zoomed in 
the inset of the figure. Its linewidth is 260 Hz in optimal 
operating conditions: A^ = A m , probe pulse of 10 ms 
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FIG. 2: Motional spectrum of the atoms in the optical lattice. 
Each point corresponds to a single measurement. Black curve: 
optimal operating conditions. The inset is a zoom on the car- 
rier. Grey curve: the motional sidebands are enhanced by 
applying long probe pulses of 200 ms at the maximum avail- 
able power at 698 nm. 

duration and 2 /iW power. The resonance plotted here 
was obtained at the maximum trapping depth and we 
observed no clear dependence of its width and contrast 
with Uq. 

The carrier is surrounded by motional sidebands 
shifted by the oscillation frequency along the lattice di- 
rection. They are hardly visible in optimal operation but 
can be enhanced by applying long pulses of 200 ms du- 
ration at the maximum available power (2 mW) of the 
probe beam. The resulting spectrum is plotted in grey 
in Fig-El The 1:5 ratio between the red and blue side- 
bands show that 80% of the atoms populate the \n z =0) 
motional state along the lattice axis, corresponding to 
a temperature of 8 /iK. The transverse temperature is 
10 /iK. 

The effect of the trap on the clock transition is mea- 
sured by locking the frequency of the probe laser to the 
carrier for various values of Uq and Al- The trap depth 
Uq is adjusted to its desired value between 200 E r and 
1400 E r by a linear ramp of 1 ms duration between the 
cooling and probing phases. This slightly decreases the 
transverse temperature while the longitudinal motion is 
adiabatically cooled by following the ramp. To lock the 
probe laser frequency we alternatively probe both sides 
of the resonance and compute an error signal from the 
difference between two successive measurements of the 
transition probability. This also gives a measurement of 
the difference between the atomic transition and the ref- 
erence cavity frequency which slowly fluctuates due to 
thermal effects. These fluctuations behave essentially as 
a sine wave of peak to peak amplitude 300 Hz and pe- 
riod 10 minutes. To derive the frequency shift due to the 
lattice, we reject the cavity frequency fluctuations by a 
factor 100 by a differential method [l5j. We interleave 
measurements at 4 different lattice depths. We run the 



3 




400 



800 1200 
U„[E] 

L r J 



N 

X 




813.40 813.42 813.44 

X [nm] 




813.32 813.36 813.40 813.44 813.48 
X [nm] 



FIG. 3: (a): Frequency shift of the clock transition vs lattice 
depth for 5 different lattice wavelength: 813.406 nm (four 
lowest points in the graph), 813.413 nm, 813.428 nm, 813.436 
nm and 813.444 nm (four highest points). Also shown is a 
linear fit of each of these five sets of data, (b): First order 
frequency shift vs Al for Uo = E r . Also plotted is a linear 
fit of these data (x 2 = 1.1). The inset shows measurements 
of A m . ■ : this work. ▲ : Ref. Q. • : Ref. 0- (c): Higher 
order frequency shift vs Al for Uo = E r . All these points 
are compatible with zero. Their average is -4(36) x 1(T 7 Hz 
(X 2 = 1-02). The arrow on the \l axis corresponds to the 
3 Po — > X Pi transition at 813.360 nm. The vertical dotted line 
is at the magic wavelength 813.428 nm. 



clock for 19 cycles before Uo is changed to the next of 
the four interleaved values. The entire sequence is re- 
peated typically 16 times. The cavity fluctuations are 
then modelled as a polynomial which is determined by a 
least square fit of the data • The data are corrected for 
the modelled cavity frequency fluctuations and avera ged . 
This yields 4 statistically independent measurements [l9| 
of the clock transition frequency versus Uq with a stan- 
dard deviation of about 5 Hz. A typical set of such points 
is shown in Fig. Of a). 

We perform a quadratic least square fit of each set 
of four points which gives a measurement of the coeffi- 
cients v\ and V2 of Eq.(Q). The coefficient v-i is plot- 
ted in Fig.^c) and found compatible with zero for the 
whole range [813.3 nm, 813.5 nm] to within a few tens of 
/zHz. A specially interesting wavelength region is around 
813.36 nm where the 3 Pq — » 1 P\ two-photon transition is 




FIG. 4: (a): Higher order frequency shift around the 3 Po — » 
3 ' F2 transition at 818.570 nm for Uo = E r . The five vertical 
arrows on the wavelength axis correspond to the hyperfine 
sub-states of 5s4/ 3 F 2 (F = 13/2 to F — 5/2 from left to 
right), (b): Atomic frequency shift vs trapping depth for two 
lattice wavelengths on both sides of the two-photon transi- 
tion to sub-state F = 13/2. The linear light shift has been 
removed for clarity. The bold line is a fit of the data with a 
parabola. 



expected. The contribution of a two-photon transition to 
V2 varies as A" 1 , with A the detuning of the lattice with 
respect to the resonance. We magnify this contribution 
by systematically spanning a frequency range of ±5 GHz 
with a 1 GHz step around the expected value (Fig.|3Jc)). 
The null results of all these measurements demonstrate 
that the higher order shift due to the 3 Po — * 1 Pi tran- 
sition is less than 1/iHz/i?;? for A^ = A m . Despite its 
proximity to the magic wavelength, this two-photon tran- 
sition is forbidden enough to not be a problem. 

This set of experiments around 813.4 nm can also 
be used to derive an accurate value of A m . Having 
shown that V2 is negligible for A^ ~ A m , better esti- 
mates of v\ are obtained with linear fits of each set of 
four points. They are plotted in Fig.|jJJb). We find 
A m = 813. 428(1) nm in agreement with previously pub- 
lished values |2j, yj as shown in the inset of FigJSfb) . The 
improvement by one order of the accuracy of this mea- 
surement is a nice illustration of the amplification of the 
effects of the lattice offered by a deep trapping potential. 

We also studied the effect of the other two photon 
transition in this wavelength region, the 3 Pq — > 3 i<2 at 
818.57 nm. When tuned 5 nm away from A m we expect, in 
addition to the effect of the two-photon coupling, a trivial 
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FIG. 5: Higher order frequency shift at the magic wavelength 
for Uo = E r . The value of the open square is the weighted 
average of the seven other points (x 2 = 1.8). 



quadratic dependence of the atomic frequency vs Uq due 
to the imperfect cancellation of v\ and to the inhomo- 
geneity of the laser intensity experienced by the atoms. 
We do observe a substantial broadening and asymmetry 
of the resonance due to this effect similar to what was re- 
ported in Ref. 16]. The associated trivial quadratic fre- 
quency shift amounts to 0.8 mHz/E^ as measured several 
GHz away on both sides of the two photon transition. 

When tuned closer to the resonance, we clearly observe 
the non-trivial quadratic frequency shift due to the two- 
photon resonance itself. The effect is shown in Fig. 01 
Quadratic frequency shifts of several kHz and changing 
sign depending on the side of the transition being probed 
are observed as shown in Fig.HJb). This is a clear sig- 
nature of the higher order effects due to the particular 
transition under investigation. When tuned exactly onto 
one of the five transitions corresponding to the hyperfine 
structure of 3 F 2 , the lattice laser induces severe loss (up 
to 90%) of atoms when the 698 nm probe laser is tuned 
to resonance. This effect, which we attribute to three- 
photon ionization from 3 Pq, was used to determine the 
position of the five hyperfine sub-states shown by arrows 
in Fig.^Ja). The hyperfine structure of 3 F 2 leads to a 
complex dependence of the quadratic frequency shift on 
the lattice wavelength around resonance. The contribu- 
tion of the five substates can interfere with each other, 
which may be the cause of the oscillations of v 2 seen in 
Fig. Ola) on both sides of the hyperfine manifold. We can 
deduce from our measurements a conservative estimate 
of the contribution of the 3 Pq — 3 F 2 resonance to the 
higher order effects at Al = A m . The grey dashed curved 
plotted in Fig.^a) scales as the inverse of the detuning 
of the lattice with respect to the center of gravity of the 
hyperfine structure of 3 _F2 and envelopes experimental 
points. When extrapolated to the magic wavelength, it 
gives a contribution to v 2 of 2 /iHz. 

Finally we have performed an extensive series of mea- 
surements of v 2 at Al = A m resulting in the data plot- 
ted in Fig.0 Their weighted average gives v 2 (\ m ) = 
7(6) /j,Hz. For a lattice depth of 10 E r the frequency shift 
is lower than 1 mHz (one sigma) or 2 x 10~ 18 in fractional 
units. This demonstrates that the frequency shift due 



to the atomic hyperpolarizability constitutes no impedi- 
ment to the accuracy of a Sr optical lattice clock down to 
the 10 -18 level. In addition, the effective laser intensity 
seen by the atoms is certainly controllable at the percent 
level |20j . The performance of the system would then be 
immune to higher order frequency shifts over a broad lat- 
tice depth range, possibly up to Uq — 100 E r . This would 
provide a powerful lever for the experimental evaluation 
at the 10~ 18 level of other effects associated for instance 
to the dynamics of the atoms in the lattice or to cold 
collisions. Collisions are expected negligible with polar- 
ized fermions, but they have to be considered if one uses 
bosonic isotopes 0, such as 88 Sr. By varying the trap- 
ping depth, one can adjust the tunnelling rate and then 
control the overlap of the wave functions of atoms con- 
fined in the lattice, allowing the study of cold collisions 
in a new regime. 
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All wavelengths throughout the paper are in vacuum. 
We use a tenth order polynomial fit of the data. For any 
polynomial order between 5 and 14, the results are un- 
changed to within a fraction of the error bars, 
except for a common offset which doesn't play any role 
in further data analysis. 

A control to within a few percents is already achieved 
in our setup as evidenced by the nice alignement of the 
measurements of v\ shown in Fig.|3fb). 



